Measurement of the Cosmic-Ray Antiproton to Proton Abundance Ratio 

between 4 and 50 GeV 
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We present a new measurement of the antiproton to proton abundance ratio, p/p, in the cosmic 
radiation. The HEAT-pbar instrument, a balloon borne magnet spectrometer with precise rigidity 
and multiple energy loss measurement capability, was flown successfully in Spring 2000, at an 
average atmospheric depth of 7.2 g/cm 2 . A total of 71 antiprotons were identified above the vertical 
geomagnetic cut-off rigidity of 4.2 GV. The highest measured proton energy was 81 GeV. We find 
that the p/p abundance ratio agrees with that expected from a purely secondary origin of antiprotons 
produced by primary protons with a standard soft energy spectrum. 

PACS numbers: 95.85.Ry, 96.40.De 



Antiprotons constitute a rare but interesting compo- 
nent of the cosmic radiation. They are secondary cosmic- 
ray particles, generated in nuclear interactions of high- 
energy cosmic rays with the interstellar medium (ISM). 
It remains an open question whether there are significant 
additional contributions that have a different and per- 
haps more exotic origin. The kinematic threshold energy 
for p production in p-p collisions of primary cosmic-ray 
protons causes a p energy spectrum, and a p/p inten- 
sity ratio, that decline rapidly from a few GeV towards 
lower energies. Solar modulation inside the heliosphcre 
softens this 'cutoff' and leads to uncertainties in the pre- 
dicted flux at low energy. Additional sources of antipro- 
tons might be evaporating primordial black holes (PBH) , 
or annihilating supersymmetric particles. The PBH con- 
tribution would be expected to be significant at energies 
well below 1 GeV, while supersymmetric particle anni- 
hilations, for instance neutralinos, could also affect the 
antiproton intensity at higher energy, above several GeV. 
At those higher energies, the p energy spectrum will be 
essentially unaffected by uncertainties due to solar mod- 
ulation. 

The energy spectrum of antiprotons measured near 
Earth carries the imprint of losses during propagation 
through the Galaxy and thus, is a sensitive probe of the 
confinement environment of protons. For instance, if the 
propagation path length A and the diffusion coefficient 
for protons depended on energy E in the same way as 
has been observed for the heavy cosmic ray nuclei (i.e. 

, one would expect the antiproton fraction p/p 
to gradually decrease with increasing energy above a few 
GeV. This behavior would be equivalent to that of the 
intensities of secondary spallation nuclei such as Li, Be, 
and B, relative to those of their heavier primary parents 
C and O. In order to obtain a self-consistent model of 



the propagation of protons in the Galaxy the observed p 
spectrum must also be compared with measurements of 
positrons and gamma rays, which also result from nuclear 
interactions in the ISM (mostly via ir + and ir° decay). 

It would be difficult to explain if the antiproton frac- 
tion p/p were found to be constant or if it increased at 
higher energy. Models of extragalactic origin are un- 
likely because of constraints in intergalactic transport 
. Closed Galaxy models have been suggested ^ which 
would boost the p/p ratio at high energy, but these would 
also predict a higher abundance of He 3 than has been ob- 
served J|, H . It has been proposed that the primary pro- 
ton spectrum in distant regions of the Galaxy is harder 
than locally observed ||. This would also lead to en- 
hanced p production in these regions, although not nec- 
essarily to an enhanced p fraction near Earth. Thus it 
appears that a constant or rising p/p fraction might in- 
deed require the presence of a primary and possibly exotic 
source of antiprotons. 

Observationally, the situation has been unclear. Mea- 
surements at low energies, in particular the series of ob- 
servations with the BESS instrument || |7), have pro- 
vided the p energy spectrum with good statistical accu- 
racy from ~ 0.2 to 3 GeV. These results are in good 
agreement with interstellar secondary production mod- 
els. PBH contributions, if they exist, are hidden by the 
uncertainties due to solar modulations. Above 5 GeV, 
the results reported in the three previous measurements 
H [| |l^, |ll| are statistics limited, and no solid conclu- 
sion about the shape of the energy dependence of the an- 
tiproton fraction can be drawn. The current HEAT-pbar 
experiment has been developed to clarify this situation 
with a series of balloon flights. 

A schematic view of the HEAT-pbar instrument is 
shown in Fig. |l|. It consists of a superconducting magnet 
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with a drift tube hodoscope at its center, combined with 
stacks of multi-wire proportional chambers above and be- 
low the hodoscope. Two layers of scintillators, one at the 
very top and one at the very bottom of the instrument, 
provide time-of-flight (ToF) measurements and, together 
with a scintillator just above the hodoscope, form the 
event trigger. 

The spectrometer has been described in detail |l2). It 
consists of 479 drift tubes in twenty four layers. It is 
mounted in the room-temperature bore of a supercon- 
ducting magnet with a 10 kG central field, and measures 
the particle trajectory within the magnetic field, provid- 
ing both the particle rigidity (momentum/charge) and 
charge sign. It has a single-point tracking accuracy of 
better than 70 /Ltm, an average track length of 58 cm, 
and typically measures 15 points along a particle trajec- 
tory. The continuous-tracking approximation then yields 
a mean Maximum Detectable Rigidity (MDR) of 170 GV, 
and the rigidities of particles up to ~ 60 GV can be re- 
liably measured. 

Antiproton flux measurements require excellent parti- 
cle identification for background discrimination. The pri- 
mary sources of background to the p flux are electrons, 
and negatively charged muons, pions and kaons produced 
in the atmosphere as well as in the material above and 
in the detector. To provide mass discrimination between 




FIG. 1: Schematic diagram of the HEAT-pbar instrument. 
The scintillator (S) above the hodoscope is part of the event 
trigger. Upper and lower ToF scintillators are 2.8 m apart. 



antiprotons and these particles we measure multiple sam- 
ples of the ionization energy loss. The logarithmic rise 
in the mean rate of energy loss for a relativistic charged 
particle is used to determine the Lorentz factor of the 
particle from which, together with the rigidity measure- 
ment, the mass is obtained. 

The multiple dE/dx detector consists of a stack of 140 
segmented multi-wire proportional chambers, each pro- 
viding a measurement of the specific ionization loss. In 
order to maximize the particle identification power, the 
proportional chambers are filled with Xenon (and 5% 
CH4), which exhibits an increase in ionization loss rate 
of 70% between minimum ionization and relativistic sat- 
uration Jl3[ . 

The HEAT-pbar instrument was flown from Ft. Sum- 
ner, NM on June 3, 2000. The detector was at float 
altitude for 22 hours, at an average atmospheric over- 
burden of 7.2 g/cm 2 . More than 1.9 million events were 
recorded over an integrated live time of 16.2 hours at a 
residual pressure between 4.5 and 8.6 mbar. The average 
vertical geomagnetic cut-off rigidity along the flight path 
is 4.2 GV. The instrument performed flawlessly during 
the flight. 

Albedo particles, which mimic antiparticles in the 
spectrometer, can be efficiently rejected with the ToF 
measurement. The flux of upward-going relativistic par- 
ticles is roughly 10~ 3 that of relativistic downward-going 
particles, and thus a rejection power of 1000 is required 
to keep the contamination of these particles below 1% in 
the final antiproton data set. The standard deviation in 
the velocity distribution for relativistic protons is 0.093c. 
This results in a rejection power against upward moving 
particles which is several orders of magnitude better than 
required. 

In analyzing events recorded by the dE/dx system, 
first suitable selections on tracking quality are made and 
events are selected for which the incident particle transits 
the entire dE/dx chamber system. The effect of high- 
energy tails in the Landau distribution is minimized by 
computing, for each event, a restricted average specific 
ionization (dE / dx res ) which is the average ionization sig- 
nal measured by 50% of the dE/dx chambers with signals 
smaller than the median for this event. Note that the 
precise value of the selected fraction is not very critical, 
but that 50% is close to optimal. Fig. 2a and b show his- 
tograms of the (dE I dx res ) response for negative and pos- 
itive particles, respectively, having rigidities in the range 
4.5-6 GV. In these distributions, the peak at smallest 
(dE I dx res ) corresponds to protons and antiprotons, the 
next peak to n + / fi + and 7r~//x~, and the peak at large 
(dE / dx res ) to e + and e~. Compared to the antiproton 
flux, the kaon production is small, although not negligi- 
ble. Our Monte Carlo simulations show that our event se- 
lection criteria reduce this contribution to negligible lev- 
els, since kaons result from interactions in or near the in- 
strument. A somewhat better mass resolution can be ob- 
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FIG. 2: Histograms of the (dE/dx rea ) response for negative 
(a) and positive (b) particles in the rigidity range 4.5 - 6 GV. 
The distributions of particle species are gaussian to better 
than four order of magnitude (as can be seen in the rising 
edge of the proton distribution for instance). The upper dis- 
tribution in each figure shows the corresponding data before 
track quality selection criteria are applied (and thus the ma- 
jority of events in the large peak at the antiproton position 
are really proton tracks out in the tails of the track parameter 
distributions.) 



in the atmosphere and in the instrument (maximum col- 
umn density 5.7 g/cm 2 above the lower ToF counter for 
a vertically traversing particle) are applied. The correc- 
tions assume that all particles that interact inside the 
gondola are rejected by our selection criteria. 

The calculated background of secondary antiprotons 
and protons produced in the atmosphere was based on 
Pfeifer et al. Q] for the antiprotons and Papini et al. |l|] 
for the protons. Interaction and annihilation losses are 
based on the measured cross sections quoted in Kuzichev 
et al. |l(J and Denisov et al. fL7| , accounting in detail 
for the total material traversed by a particle in passing 
through the atmosphere, aluminum shell, and detector 
material. The number of antiprotons and protons in each 
energy bin obtained after applying all of these corrections 
are shown in Table |, along with the resulting p/p ratios. 
The errors quoted in this table are purely statistical. Sys- 
tematic errors resulting from uncertainties in correcting 
the particle numbers in the instrument to the top of the 
atmosphere and in the background due to particle mis- 
identification are estimated to be less than 4% of the p/p 
ratio. 

Our results are shown in Fig. along with previous 
measurements by others, and a number of recent cal- 
culations of the p/p ratio resulting from secondary pro- 
duction of antiprotons in the interstellar medium. Only 
recent measurements have been included in this figure 
§0,1,1,0,0. 

Many predictions for the p/p ratio have been published 
over the years. We show here theoretical curves that are 
consistent with the now well measured flux ratio in the 
low energy region around 1 GeV. The result of calcu- 
lations by Simon et al. |l8| are shown in the figure as 
a shaded band. The calculations are based on the leaky 
box model and the uncertainties in the flux prediction, re- 
flected by the band in the figure, are primarily uncertain- 
ties in the galactic path length distribution. The dashed 



tained by properly accounting for the dependence of the 
(dE / dx res ) signal on the rigidity within a given rigidity 
interval, but these histograms, which are representative 
of the equivalent histograms at higher energy, demon- 
strate clearly that particle identification is achieved for 
both antiprotons and positrons. The large sampling for 
each energy loss measurement produces highly gaussian 
distributions, and thus for each energy interval, the p/p 
ratio can be obtained from fits to the restricted average 
dE/dx distributions, such as those shown in Fig. 2a and 
b for the rigidity interval from 4.5 - 6 GV. 

The results are summarized in Table |. In order to 
obtain the number of protons and antiprotons at the top 
of the atmosphere, we correct for particle production in 
the atmosphere (total average column density 7.2 g/cm 2 ) 
above the instrument. In addition, corrections for inter- 
action and annihilation losses of protons and antiprotons 



TABLE I: Event selection results and p/p ratios (in 10 -4 ). 
R is the measured rigidity at the spectrometer and T is the 
corresponding kinetic particle energy at the top of the at- 
mosphere. N p and Np are the number of observed protons 
and antiprotons for each energy bin, respectively. Np orr ' and 
N^-° rr ' are the extrapolated number of protons and antiprotons 
at the top of the atmosphere. The pion/muon background 
due to tails in the (dE/dx res ) distributions in the five rigidity 
bands is (0.2, 0.3, 0.4, 0.7, 0.8) counts and is included in the 
corrections. 



R (GV) T (GeV) 



N, 



% Np° rr - Nf r 



ratio 



4.5 - 6.0 3.7 
6.0-10.0 5.1 
10.0 - 15.0 9.1 - 
15.0 - 25.0 14.1 
25.0 - 50.0 24.1 



5.1 119361 18 124814 

9.1 141447 23 148952 

14.1 60727 21 64971 

24.1 37742 15 40141 

49.1 8773 1 9090 



0.50 
0.39 
+0.46 
0.37 
+ 1.01 
0.81 

3.21+ 142 



13.9 1.11 
16.9 1.13 
18.9 2.91 
12.9 10 
< 2.1(90%) 
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FIG. 3: Compilation of observed p/p flux ratios at the top of 
the atmosphere, compared with model calculations for sec- 
ondary andprimary antiproton production: BESS 95&97 
§, BESS §,IMAX §, MASS91 §, CAPRICE94 jto|, 
CAPRICE98 Q. The calculations of the p/p ratio are from 
H (MSR-1, MSR-2) and Jl|](SMR). Possible primary contri- 
butions to the p/p spectrum arising from evaporating primor- 
dial black holes ]1S| (MMO) and from neutralino annihilation 
[po[ (J&K) are also shown. 



and solid line in Fig 
Moskalcnko et al. IB 



'J show the results of calculations by 
within a self-consistent CR propa- 
gation model. The dashed line represents the case of a 
proton injection spectrum that is much harder than lo- 
cally observed, which has been proposed to explain the 
observed high continuum gamma-ray emission above ~ 
1 GeV Q. A standard proton injection spectrum, con- 
sistent with the locally observed one, is reflected in the 
solid line. The sensitivity of the p/p ratio to the nucleon 
injection spectrum above a few GeV makes antiproton 
measurements at energies above a few GeV an important 
test for CR models. Our data are in good agreement with 
the 'standard spectrum' calculations || at high energy, 
and do not support an antiproton to proton ratio ap- 
proaching 10 -3 at energies above 20 GeV, in contrast to 
recent CAPRICE measurements [[n] . Our result does not 
support models which are based on hard nucleon injec- 



tion spectra. At energies covered by the measurements 
presented here, secondary p production with a nucleon 
injection spectrum consistent with the locally observed 
one describes the data well. 

The HEAT-pbar instrument is scheduled for additional 
balloon flights and we expect to statistically improve the 
data and to further clarify the experimental situation. 
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